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RESEARCH  SUMMARY 

Fine  fuel  samples  from  Northwestern  forest  areas  were 
subjected  to  an  electric  arc  discharge  that  simulated  the 
lightning  continuing  current.  The  important  fuel  state  para- 
meters and  the  ignition  probabilities  were  determined  using 
logistic  regression. 

The  ignition  probabilities  for  duff  from  short-needled  coni- 
fer species  were  found  to  depend  almost  entirely  on  the 
depth  of  the  fuel  bed.  In  the  case  of  litter  and  duff  from  long- 
needled  species,  ignition  probabilities  depend  mostly  on  the 
ovendry  fuel  moisture. 

Ignition  probabilities  also  depend  on  the  duration  of  the  arc 
only,  with  current  flow  and  diameter  having  no  detectable 
effect.  Because  the  continuing  current  durations  can  be  fit 
well  by  a  Weibull  distribution  for  both  positive  and  negative 
discharges  to  ground,  the  probabilities  conditional  on  the 
current  duration  can  be  converted  to  marginal  probabilities 
per  discharge.  The  necessary  integrations  were  done  and 
the  results  are  presented. 

Marginal  ignition  probabilities  are  given  in  four  forms. 
Tables  give  the  probabilities  for  selected  values  of  the  inde- 
pendent variables.  Graphs  of  the  probabilities  are  also 
shown.  Equations  based  on  curve  fits  to  the  tabulated  val- 
ues accompany  the  graphs.  Finally,  a  FORTRAN  program 
that  uses  a  rapid  integration  scheme  to  give  the  probabilities 
is  included  in  subroutine  form.  All  these  products  are  valid 
only  over  the  range  of  values  for  the  independent  variables 
as  specified  in  the  paper. 
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INTRODUCTION 

Land  management  agencies  are  now  using  lightning 
location  equipment  in  fire  management  operations.  Some 
implications  of  this  technology  were  set  forth  by  Latham 
(1979),  including  the  steps  necessary  to  accomplish  the 
location  and  timing  of  lightning-caused  fires.  One  of  these 
steps  is  assessing  the  probability  that  lightning  discharge 
to  ground  will  ignite  fuels. 

A  previous  investigation  of  ignition  probability  (Fuquay 
and  others  1979)  used  an  energy  argument  in  which  the 
available  energy  density  in  the  lightning  continuing  cur- 
rent was  compared  to  the  fuel's  heat  of  ignition  as  used  in 
fire  models.  The  authors  concluded  that  further  investiga- 
tion was  needed  to  improve  the  accuracy  of  the  ignition 
probability  algorithm. 

The  details  of  that  further  investigation  are  contained 
in  a  series  of  papers.  The  first  paper  (Latham  1987)  de- 
scribes the  equipment  needed  to  do  the  testing.  The  sec- 
ond (Latham  1989)  reports  the  details  of  the  electric  arc 
as  an  energy  source  as  it  was  used  in  the  experiments. 
The  third  and  present  publication  describes  the  ignition 
experiments  and  the  results.  This  publication  also  for- 
mats the  results  for  application  and  discusses  probability 
of  ignition  for  lightning  continuing  currents.  The  report 
provides  tables,  curves,  formulas,  and  algorithms  for  im- 
plementation of  the  ignition  probabilities  together  with 
suggestions  for  use. 

DATA  AND  ANALYSIS 

The  electric  arc  source  was  used  to  penetrate  1,666  fuel 
beds  with  a  dc  arc  discharge.  The  current  and  duration  of 
the  discharge  were  varied,  as  were  the  fuel  bed  composi- 
tion, depth,  packing  ratio,  and  moisture.  For  each  fuel 
bed,  the  result  of  ignition  or  no  ignition  was  recorded.  The 
fuel  types  used  in  the  analysis  were  ponderosa  pine  litter, 
punky  wood  (both  in  chunks  and  fine  powder),  peat  moss 
(commercial  samples),  Douglas-fir  duff,  lodgepole  pine 
duff,  and  high-altitude  mixed  species  (primarily  Engel- 
mann  spruce)  duff. 

The  dependent  variable  has  only  two  values  (1  =  ignition 
or  0  =  no  ignition),  so  a  logistic  regression  model  was  used 
to  estimate  the  probability  of  ignition.  The  logistic  model 
gives  a  number  between  0  and  1,  which  is  appropriate 
since  we  are  estimating  a  probability.  The  maximum 
likelihood  estimates  for  model  parameters  were  obtained 


fi-om  the  SAS  procedure  CATMOD.  The  models  were 
constructed  fi-om  the  original  observations  rather  than  by 
grouping  the  data. 

One  step  in  the  investigation  was  determining  which 
independent  variables  to  include  in  the  models.  The  sta- 
tistic used  for  these  tests  was  minus  twice  the  maximized 
log  likelihood  (Efron  1978).  When  candidate  variables  are 
added  one  at  a  time,  the  difference  in  log  likelihoods  has  a 
chi-square  distribution  with  1  degree  of  freedom. 

After  the  models  were  estimated,  a  chi-square  type 
statistic,  Cg,  was  used  to  evaluate  how  well  the  models  fit 
(Lemeshow  and  Hosmer  1982).  The  statistic  is  calculated 
from  a  table  of  observed  and  expected  values  and  has  a 
distribution  closely  approximated  by  a  chi-square  distri- 
bution with  8  degrees  of  freedom.  Test  results  by  fuel 
type  are  shown  in  table  1.  A  p-value  greater  than  0.10 
indicates  that  the  model  fits  the  data  well. 

Two  important  results  were  apparent  from  the  analy- 
ses. The  primary  result  was  that  the  most  important  of 
the  variables  associated  with  the  discharge  was  its  dura- 
tion. We  checked  this  carefully,  because  of  our  original 
conjecture  that  the  energy  density  would  be  important. 


Table  1 — Evaluation  of  the  fit  of  logistic  regression  equations 
to  the  ignition  data.  Cg  and  p-value  are  explained  in 
the  text 


Fuel 


Sample  size  Cg 


p-value 


Ponderosa  pine 

litter  722  8.61  0.38 
Punky  wood 

(rotten,  chunky)  136  4.01  .86 
Punky  wood  powder 

(4.8  cm  deep)  105  4.43  .82 
Punky  wood  powder 

(2.4  cm  deep)  70  2.75  .95 
Lodgepole  pine 

(duff)  90  6.81  .56 
Douglas-fir 

(duff)  117  1.62  .99 
Engelmann  spruce 

(duff)  160  11.86  .16 
Peat  moss 

(commercial)  266             3.11  .93 
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The  fact  that  the  duration  alone  can  characterize  the 
discharge  for  ignition  is  very  significant.  It  is  the  one 
variable  associated  with  the  continuing  current  discharge 
that  is  relatively  easy  to  measure  in  the  field  (as  opposed 
to  the  current  or  the  diameter).   Second,  we  found  that 
the  moisture  content  (for  moistures  less  than  40  percent) 
plays  a  very  small  role  in  the  ignition  of  duff  fi^om  short- 
needled  species;  the  depth  of  the  bed  is  more  influential. 

LOGISTIC  REGRESSION 
COEFFICIENTS 

We  broke  the  coefficients  of  the  exponential  in  the 
logistic  regression  for  each  fuel  type  into  two  parts  for 
convenience  in  later  processing.  One  coefficient  (B)  multi- 
plies the  duration  of  the  arc  discharge,  and  the  other 
(A)  contains  the  coefficients  and  variables  that  are  specific 
to  the  fuel,  such  as  fuel  moisture,  bed  depth,  bulk  density, 
and  so  on.  These  coefficients  are  used  to  determine  the 
probability  of  ignition  using  the  equation: 

pc\(A,B,t)  =  1/(1  +  exp(-A  -Bt))  (1) 

where  the  experimentally  determined  coefficients  are 
given  in  table  2.  These  values  are  used  with  equation  (1) 
to  generate  conditional  ignition  probabilities  when  the 
duration  of  a  continuing  current  event  and  fuel  parame- 
ters are  known. 

Present  lightning  location  systems  do  not  measure  the 
duration  of  the  continuing  current.  In  fact,  they  do  not 
even  identify  those  lightning  discharges  that  have  a  con- 
tinuing current  component.  This  situation  can  be  partly 
corrected.  We  do  have  a  number  of  measurements,  at 
least  for  northwestern  storms,  for  the  duration  of  continu- 
ing currents  from  both  negative  and  positive  discharges  to 
ground.  These  measurements  can  be  combined  with  the 
conditional  probabilities  to  give  a  probability  per  light- 
ning event. 


MARGINAL  IGNITION  PROBABILITY 

What  we  want  to  know  is  the  probabihty  of  an  ignition, 
given  the  occurrence  of  an  electrical  discharge  in  the  fuel, 
that  is,  the  marginal  probability  of  ignition  with  respect  to 
discharge  duration.  We  can  now  estimate  conditional 
ignition  probability,  given  the  duration  of  the  discharge 
and  the  fuel  state,  ft-om  equation  (1)  and  table  2.  If  we 
know  the  probability  distribution  for  continuing  current 
duration,  we  can  find  the  marginal  probability  of  ignition 
by  multiplying  the  two  probabilities  and  adding  them  up 
over  all  durations. 

We  have  determined  that  the  continuing  current 
probability  is  best  fit  by  the  WeibuU  (see  appendix  A) 
distribution  (Kappenman  1988): 

p(«.n,s)  =  (J)(^)""'  exp[-(Vs)"]  (2) 

where  t  is  the  duration  (msec)  and  the  values  of  n  and  s 
are  1.6  and  207.8  for  negative  discharges,  and  2.3  and 
69.3  for  positive  discharges. 

The  adding-up  is  done  by  the  integral: 

oo 

pi(n.s,A,B)  =  ^p(t,n,s)  pci  (A,B,t)dt  (3) 
0 

that  must  be  solved  by  a  numerical  integration.  We  per- 
formed the  integrations  necessary  and  obtained  marginal 
probability  distributions  as  given  in  table  3  and  appendix 
B.  The  simple  equations  of  table  3  result  from  a  "best  fit" 
approach  starting  with  simple  functions. 

We  have  also  developed  a  FORTRAN  algorithm  (appen- 
dix C)  for  the  marginal  probability.  This  algorithm  per- 
forms the  integration  necessary  to  obtain  the  marginal 
distribution  as  presented  above.  The  integration  is  a  6- 
point  Gaussian  quadrature  using  Laguerre  polynomials. 
The  results  check  within  2  to  3  percent  of  the  Romberg 
quadrature  used  to  calculate  the  values  of  appendix  B. 
This  precision,  considering  the  other  uncertainties  in 
application  of  the  algorithms,  is  entirely  adequate. 


Table  2 — Coefficients  for  tfie  probability  distributions  of  Equation  (2).  Mf  is  the  fuel 

moisture  in  percentage  of  dry  weight,  5  is  the  duff  depth  in  cm,  and  pfa  is  the 


bulk  density  in  g/cm^ 

Fuel                                               A  B 

Ponderosa  pine  litter  0.97-0. 19M  0.012 

Punky  wood  (rotten,  chunky)  -0.59  -  OASMf  .005 

Punky  wood  powder  (4,8  cm  deep)  1 .2  -  0.1 2M/  .002 

Punky  wood  powder  (2.4  cm  deep)  0.13  -  O.OSMf  .005 

Lodgepole  pine  (duff)  -5.6  +  0.685  007 

Douglas-fir  (duff)  -7.1  +  1.45  006 

Engelmann  spruce  (duff)  0.79  -  0.081  Mf- 8. 5piJ  .011 

Peat  moss  (commercial)  0.42- 0.1 2Mf  .005 
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Table  3 — Equations  approximating  the  curves  of  figs.  3  through  10.  Mis  the  fuel  moisture  in  percentage  of 
dry  weight  and  5 is  the  bed  depth  in  cm.  These  equations  are  valid  only  within  the  specified 
ranges 


Fuel 


pineg 


pipes 


Ponderosa  pine  litter 
Punky  wood  (rotten,  chunky) 
Punky  wood  powder  (4.8  cm  deep) 
Punky  wood  powder  (2.4  cm  deep) 
Lodgepole  pine  (duff) 
Douglas-fir  (duff) 
Engelmann  spruce  (duff) 
Peat  moss  (commercial) 


1.04  exp(-0.054MO 
0.59  exp(-0.094/lf/) 
0.9  exp(-0.056Aff) 
0.73  -0.011  M 
(1+exp(3.84-0.65))-' 
(1+exp(5.48-1.286))- 
0.8- 0.01 4/lf/ 
0.84exp(-0.06Aff) 


0.92  exp(-0.087/W/) 
0.44  exp(-0.11A^ 
0.86  exp(-0.06AW) 
0.6-0.11Mf 
(1+exp(5.13-0.685))-' 
(1+exp(6.69-1.398))-' 
0.62  exp(-0.05M/) 
0.71  exp{-0.07Mf) 


Note:  Range  of  validity  is  0  to  40  percent  for  Mf  and  0  to  10  cm  for  5.  Do  not  use  outside  this  range. 


SUGGESTIONS  FOR  USE 

Ignildon  probabilities  have  a  wide  range  of  uses,  from 
preseason  planning  to  dispatch  of  fire  suppression  forces. 
Specific  applications  are:  daily  fire-danger  rating, 
preattack  (short-range)  planning,  prepositioning  of  fire- 
fighting  forces,  and  dispatch  of  fire-spotting  flights. 
These  uses  of  ignition  probability  depend  on  the  available 
knowledge  about  the  variables  involved  in  the  calcula- 
tions. Generally,  the  farther  from  real-time  the  applica- 
tion is,  the  poorer  the  knowledge  about  the  fuel  variables. 

For  all  applications  that  are  not  real-time,  the  marginal 
probability  distributions  should  be  used  with  area  light- 
ning statistics.  If  the  duration  of  discharge  is  not  known, 
but  the  type  of  discharge  is  (as  in  a  real-time  application 
using  lightning  location  systems),  the  marginal  probabili- 
ties given  in  appendixes  B  and  C  should  be  applied.  The 
simple  formulas  in  table  3  may  be  used  for  this  purpose  as 
well. 

When  applying  tabulated,  graphed,  or  simple  formula 
approximations,  the  probabilities  given  are  per  continuing 
current.  Almost  all  positive  discharges  have  a  continuing 
current.  The  probabilities  for  positive  discharges  thus  do 
not  need  to  be  corrected,  and  may  be  used  as  is  per  posi- 
tive event.  Negative  discharges  to  ground  do  not  all  have 
continuing  currents.  In  fact,  only  about  one  in  five  do. 
When  using  the  probabilities  pneg,  then,  a  correction 
must  be  applied;  that  is,  a  multiplier  of  0.2.  The 
FORTRAN  program  of  appendix  C  has  this  correction 
built  in,  so  if  the  program  is  checked  against  the  tables  (as 
it  should  be  when  implemented),  the  negative  probabili- 
ties will  seem  to  be  too  small. 


Application  of  the  distributions  to  species  other  than 
those  that  we  studied  is  possible  according  to  the  sugges- 
tions of  table  4.  These  are  based  on  Agriculture  Hand- 
book 519  (Little  1978)  and  consultation  with  Dr.  J.  K. 
Brown  of  the  Intermountain  Fire  Sciences  Laboratory. 
In  addition,  we  feel  that  short  grasses  will  behave  much 
the  same  as  the  loose,  low-organic-content  peat  moss  that 
we  used  in  the  experiments,  and  tall  grasses,  such  as 
sawgrass,  will  behave  as  ponderosa  pine  litter,  because 
the  ignition  will  take  place  in  the  clumps  of  dead  grass 
rather  than  in  the  standing  dead  spears.  If  species  are 
not  mentioned  in  table  4,  we  suggest  that  consideration  of 
cover  type  be  the  guide. 

Because  of  species/cover  type  dependence,  the  primary 
use  of  probabilities  will  eventually  be  in  a  geographic 
information  system  (GIS).  In  such  use,  a  calculation  for 
each  discharge  may  not  be  necessary.  Marginal  ignition 
probabilities  can  be  calculated  in  advance  of  storms.  To 
do  this,  of  course,  one  needs  the  relative  number  of  posi- 
tive and  negative  strikes.  (About  one  positive  for  about 
every  50  negatives.)  Indeed,  for  several  species,  the  igni- 
tion probability  is  nearly  static,  because  these  depend 
primarily  on  duff  depth.  The  ignition  probabilities  can  be 
displayed,  perhaps  in  color,  and  the  discharges  superim- 
posed on  them.  The  resulting  display  would  be  projected 
ignitions,  not  surviving  fires.  In  order  to  determine  sur- 
vivability, rain  and  projected  fire  growth  models  should  be 
connected  to  ignitions.  A  learning  program  operating  in 
the  GIS  should  be  possible  to  improve  the  usefulness  of 
the  ignition  probabilities  and  survival  potential.  More 
specific  application  of  ignition  probabilities  depends  on 
identification  of  continuing  currents  and  measurement 
of  their  duration  in  lightning  locating  systems. 
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Table  4 — Application  of  fuels  studied  to  other  species 


Studied  fuels 


Ponderosa  pine 

Lodgepole  pine 

Douglas-fir 

Engelmann  spruce 

Peat  (commercial) 

Longleaf  pine 

Eastern  white  pine 

Bald  cypress 

Spruce- fir 

Ail  grasses 

Slash  pine 

Limber  pine 

All  spruce  (except  Engelmann) 

All  hemlock 

Eastern  tamarack  (peat) 

Loblolly  pine 

Tamarack 

Giant  sequoia 

Subalpine  fir 

Black  spruce  (peat) 

Pond  pine 

Jack  pine 

Eastern  hemlock 

Pacific  silver  fir 

Other  peat 

Pitch  pine 

Virginia  pine 

All  cedars  and  juniper 

Grand  fir 

All  mosses 

Shortleaf  pine 

Sand  pine 

Noble  fir 

Red  pine 

Spruce  pine 

California  red  fir 

K^nntorov/  nino 
iviui  1  ic^i  c^y  ic 

Pinunnc 
1  II  ly  ui  1  o 

Sugar  pine 

White  fir 

Jeffrey  pine 

Redwood 

Digger  pine 

Western  white  pine 

Knobcone  pine 

Bristlecone  pine 
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APPENDIX  A:  WEIBULL  DISTRIBUTIONS  FOR  LIGHTNING  CURRENTS 


Continuing  current  duration  distributions  are  well  fit  by  the  Weibull 
distribution: 


p(t,n,s)  := 


M] 


n-l  - 
•e 


(1) 


where  t  is  time,  n  is  the  shape  parameter,  and  s  the  scale  parameter. 

The  Weibull  distribution  is  the  statistic  of  choice  for  systems  in  which 
failure  of  the  whole  is  contingent  on  failure  of  a  part  or  peirts  of  the 
whole.  The  probability  of  failure  of  the  whole  increases  with  time  accord- 
ing to  the  value  of  n.  Previously,  the  distribution  of  continuing  currents 
was  fit  by  a  lognormal  distribution.  We  found  the  Weibull  to  be  a  more 
satisfactory  fit. 

The  distributions  for  negative  and  positive  lightning  strokes  to  ground 
(beised  on  data  fi-om  141  negative  strokes  and  54  positive  strokes)  are 
found  fi-om  equation  (1)  with  the  values: 


pneg  :  =  p 


^  T,1.6,207.8j 


ppos  :  =  p 


T,2.3,  69.3 


(2) 


where  x  is  duration  in  msec,  pneg  is  the  probability  of  occurrence  of  a 
negative  continuing  current  of  duration  -j,  and  ppos  is  the  probability  of 
occurrence  of  a  positive  continuing  current  of  duration  x-  The  results  of 
this  exercise  are  given  in  table  5  and  figures  1  and  2. 


Table  5 — Values  of  the  Weibull  distribution  for 
continuing  currents.  Pneg  for 
negative  discharges  to  ground,  and 
ppos  for  positive  discharges  to 
ground,  x  is  the  continuing  current 
duration  in  msec 


Weibull  distributions 

X 

Dnea(XlOOO) 

DDOS(XIOO) 

0 

0 

0 

10 

1.237 

0.265 

20 

1.846 

.623 

30 

2.304 

.966 

40 

2.667 

1.225 

50 

2.957 

1.354 

60 

3.186 

1.342 

70 

3.363 

1.208 

80 

3.495 

.995 

90 

3.586 

.752 

100 

3.640 

.523 

110 

3.662 

.335 

120 

3.656 

.198 

130 

3.624 

.107 

140 

3.570 

.054 

150 

3.497 

.025 

160 

3.408 

.010 

170 

3.305 

.004 

180 

3.191 

.001 

190 

3.068 

200 

2.938 

210 

2.802 

220 

2.664 

230 

2.524 

240 

2.383 

250 

2.243 

260 

2.105 

270 

1.970 

280 

1.838 

290 

1.710 

300 

1.587 

310 

1.469 

320 

1.357 

330 

1.249 

340 

1.148 

350 

1.052 

360 

.963 

370 

.878 

380 

.800 

390 

.727 

400 

.659 

410 

.596 

420 

.538 

430 

.485 

440 

.436 

450 

.391 

460 

.351 

470 

.314 

480 

.280 

490 

.249 

5 


.005  \ 


Ox  500 

Figure  1 — A  plot  of  the  Weibull  probability  distribution 
(pneg)  for  negative  discharges  to  ground,     is  the  dura- 
tion in  msec. 
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APPENDIX  B:  GRAPHS,  TABLES,  AND  SIMPLE  EQUATIONS  FOR  MARGINAL 
IGNITION  PROBABILITIES 


pipos  =  0.92  exp(-0.087Mf) 


Figure  3 — Marginal  ignition  probabilities  for  ponderosa  pine  litter. 
Pineg  is  the  probability  of  ignition  for  negative  discharges,  pipos  for 
positive  discharges.  Tabular  values  and  simplified  equations  are 
given  for  convenience. 
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pipos  =  0.439  exp(-0.114Mf) 

Figure  4 — Marginal  ignition  probabilities  for  rotten,  chunky  punky 
wood.  Pineg  is  the  probability  of  ignition  for  negative  discharges, 
pipos  for  positive  discharges.  Tabular  values  and  simplified  equa- 
tions are  given  for  convenience. 
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pineg  =  0.9  exp(-0.056Mf) 


pipos 
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40 


pipos  =  0.86  exp(-0.061Mf) 


Figure  5 — Marginal  ignition  probabilities  for  powdered  punky  wood, 
depth  4.8  cm.  Pineg  is  the  probability  of  ignition  for  negative  dis- 
charges, pipos  for  positive  discharges.  Tabular  values  and  simpli- 
fied equations  are  given  for  convenience. 
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pineg  =  0.73  -0.011  Mf 


0  Mf  40 


pipos  =  0.6-0.11Mf 

Figure  6 — Marginal  ignition  probabilities  for  powdered  punky  wood, 
2.4  cm  deep.  Pineg  is  the  probability  of  ignition  for  negative  dis- 
cfiarges,  pipos  for  positive  discharges.  Tabular  values  and  simpli- 
fied equations  are  given  for  convenience. 
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pipos  =  1/(1  +  exp(5.13-0.685)) 

Figure  7 — Marginal  ignition  probabilities  for  lodgepole  pine  duff. 
Pineg  is  tfie  probability  of  ignition  for  negative  discharges,  pipos  for 
positive  discfiarges.  Tabular  values  and  simplified  equations  are 
given  for  convenience. 
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pipos  =  1/(1  +  exp(6.69-1.395)) 

Figure  8 — Marginal  ignition  probabilities  for  Douglas-fir  duff.  Pineg 
is  thie  probability  of  ignition  for  negative  discharges,  pipos  for  posi- 
tive discharges.  Tabular  values  and  simplified  equations  are  given 
for  convenience. 
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pineg  =  0.8  -0.014Mf 


pipos  =  0.62  exp(-0.051Mf) 


Figure  9 — Marginal  ignition  probabilities  for  high-altitude  mixed  species, 
primarily  Englemann  spruce.  Pineg  is  the  probability  of  ignition  for  nega- 
tive discharges,  pipos  for  positive  discharges.  Tabular  values  and  simpli- 
fied equations  are  given  for  convenience. 
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pipos  =  0.715  exp(-0.072Mf) 

Figure  10 — Marginal  ignition  probabilities  for  commercial  peat 
moss.  Pineg  is  the  probability  of  ignition  for  negative  discharges, 
pipos  for  positive  discharges.  Tabular  values  and  simplified  equa- 
tions are  given  for  convenience. 
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APPENDIX  C:  A  FORTRAN  PROGRAM  FOR  ELECTRIC  ARC  IGNITION 


c 

SUBROUTINE  TH0R5  ( I FUEL, FM , DEPTH , I STROKE , PROS) 

C 

C  *********************************************** 

c 

C  TO  CALCULATE  THE  MARGINAL  IGNITION  PROBABILITY  GIVEN  A  NATURAL  ELECTRICAL 
C  DISCHARGE  OF  TYPE  <STROKE>  INTO  FUEL  TYPE  <FUEL>  WITH  MOISTURE 
C  <FM>  %  AND/OR  DEPTH  <DEPTH>  IN  CM.   NOTE  THAT  NOT  ALL  NEGATIVE  DISCHARGES 
C  HAVE  A  CONTINUING  CLTIRENT !  !  ! 

C  THE  SUBROUTINE  USES  GAUSSIAN  QUADRATURE,  WITH  THE  COEFFICIENTS  FOR 
C  THE  LAGUERRE  POLTOOMIALS  FROM: 


C  "HANDBOOK  OF  MATHEMATICAL  FUNCTIONS" 

C  M.  ABRAMOWITZ  AND  I.   STEGUN,  EDS. 

C  US  GOVT.  PRINTING  OFFICE,  APS -55  1964 

C 


C  THIS  PROGRAM  IS  A  PART  OF  "IGNITION  PROBABILITIES  FOR  WILDLAND  FUELS" 

C    BY  DON  LATHAM  AND  JOYCE  SCHLIETER. 

C 

C  !  !  ! ! ! !D0  NOT  ALTER  ANY  NUMBERS  IN  THIS  SUBROUTINE! ! ! ! ! ! 
C 

C  *************************************************************** 
C 

C  THE  INPUT  VARIABLE  <FUEL>  HAS  THE  VALUE  1-8  ACCORDING  TO: 


C  1)  PONDEROSA  PINE  LITTER, 

C  2)  ROTTEN,  CHUNKY,  PL^NKY  WOOD, 

C  3)  PUNKY  WOOD  POWDER,  DEEP  (4.8  CM) 

C  4)  PUNKY  WOOD  POWDER,   SHALLOW  (2.4  CM) 

C  5)  LODGEPOLE  PINE  DUFF 

C  6)  DOUGU^S  FIR  DUFF 

C  7)  HIGH  ALTITUDE,  MIXED,  MAINLY  ENGLEMANN  SPRUCE 

C  8)  PEAT  MOSS  (COMMERCIAL) 

C 


C  OTHER  FUELS  WILL  USE  ONE  OF  THESE  AS  SURROGATE;  CONSULT  THE 

C  ACCOMPANYING  TEXT 

C 

C  THE  INPUT  VARIABLE  <FM>  IS  THE  FUEL  MOISTURE  IN  % 

C  THE  INPUT  VARIABLE  <DEPTH>  IS  THE  DUFF  DEPTH  IN  CM. 

C  THE  INPUT  VARIABLE  <STROKE>  IS  +1  FOR  POSITIVE  DISCHARGES, 

C  -1  FOR  NEGATIVE  DISCHARGES. 

C  THE  OUTPUT  IS  <PROB>,  THE  PROBABILITY  OF  IGNITION  PER  DISCHARGE. 

C  AGAIN,  NOTE  THAT  ABOUT  1  IN  5  NEGATIVE  DISCHARGES  HAS  A  CONTINTIING 

C  CURRENT.  THIS  IS  ACCOUNTED  FOR  IN  THIS  SUBROUTINE. 
C 

C  *************************************************************** 

c 

C  SETUP  AND  PARAMETERS 
C 

DOUBLE  PRECISION  WT ,XI , ACl , AC2 , AC3 , AC4 , B , SUM , T , Y , FX , S 

C 

C  EVERYTHING  DOUBLE  PRECISION  INTERNALLY 
C 

DIMENSION  WT(6) ,XI(6) ,AC1(8) .AC2(8) ,AC3(8) ,AC4(8) ,B(8) 

C 
C 

C 

C  WEIGHTS  FOR  ORDER  6  POLYNOMIAL 
C 

DATA  WT/0. 458964673950, 0.417000830772, 0.113373382074, 
+  1.03991974531E-2,2.61017202815E-4,8.98547906430E-7/ 

C 
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C  ZEROES  OF  ORDER  6  POLYNOMIAL 
C 

DATA  XI/0. 222846604179, 1.188932101673, 2. 992736326059, 
+  5.775143569105,9.837467418383,15.982873980602/ 

C 

C  EVERY  POSITIVE  DISCHARGE  HAS  A  CONTINUING  CURRENT,   ONLY  1  IN  5 

C  NEGATIVE  DISCHARGES  DOES. 

C 

DATA  CCPERPOS.CCPERNEG  /I., 0.2/ 

C 

C  ONLY  TYPE  7  USES  BULK  DENSITY;  WE  USE  THE  AVERAGE  FOR  THAT  FUEL 
C 

DATA  RHOB  /0.127/ 

C 

C  COEFFICIENTS  FOR  THE  FUEL  TYPES 
C 

DATA  AC1,AC2,AC3,AC4,B 
+/         0.97,   -0.59,  1.2,  0.13,   -5.6,  -7.1,  0.79,  0.42, 
+        -0.19,   -0.15,   -0.12,   -0.05,  0.0,  0.0,   -0.081,  -0.12, 
+        0.0,  0.0,  0.0,  0.0,  0.68,  1.4,  0.0,  0.0, 
+        0.0,  0.0,  0.0,  0.0,  0.0,  0.0,   -8.5,  0.0, 
+        0.012,  0.005,  0.002,  0.005,  0.007,  0.006,  0.011,  0.005/ 

C 
C 

C  THESE  ARE  THE  COEFFICIENTS  IN  THE  WEIBULL  DISTRIBUTION  FOR  CCS 
C 

DATA  ETANEG,ETAPOS,SNEG,SPOS  /I . 6 , 2 . 3 , 207 . 8 , 69 . 3/ 

C 

C  DETERMINE  DISCHARGE  TYPE 
C 

IF  (ISTROKE.EQ.l)  THEN 
ETAM-1 . /ETAPOS 
S-SPOS 

COEF-CCPERPOS 
ENDIF 

IF  (ISTROKE.EQ. -1)  THEN 
ETAM-1 . /ETANEG 
S-SNEG 

COEF-CCPERNEG 
ENDIF 

C 

C  CALCULATE  F(X)  AT  XI  AND  INTEGRATE 
C 

SUM=0 . 0 

DO  100  1-1,6 

T-S*(XI(I)**ETAM) 

Y-ACl ( IFUEL) +AC2 ( IFUEL) *FM+AC3 ( IFUEL) *DEPTH+AC4 ( IFUEL) *RH0B 
+  +B(IFUEL)*T 
FX-1./(1.+DEXP(-Y)) 
100  SUM=SUM+WT(I)*FX 
PR0B=SUM*C0EF 

C 

c 

RETURN 
END 

C 

c  ********************************************** 

c 
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Ignition  of  wildland  fine  fuels  by  lightning  was  simulated  with  an  electric  arc  discharge 
in  the  laboratory.  The  results  showed  that  fuel  parameters  such  as  depth,  moisture  con- 
tent, bulk  density,  and  mineral  content  can  be  combined  with  the  duration  of  the  simu- 
lated continuing  current  to  give  ignition  probabilities.  The  fuel  state  parameters  of  impor- 
tance and  the  ignition  probabilities  were  determined  using  logistic  regression.  Graphs, 
tables,  formulas,  and  a  FORTRAN  computer  program  are  given  for  field  use. 
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